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In this work we semiclassically analyzed the high lying eigenstates of a mixed type Hamil- 
tonian system. For the regular states we employ the Einstein-Brillouin-Keller quantization, 
while for the chaotic states, following the principle of uniform semiclassical condensation, we 
obtain a prediction for their wavefunction autocorrelation function. 

§1. Introduction 

Classical Hamiltonian systems can range from completely integrable ones, such 
as the Kepler problem and the harmonic oscillator, to fully chaotic systems (e.g. 
Bunimovich stadium, Sinai billiards). Yet the majority of classical systems fall into 
neither category but are of the generic, mixed type. The phase space of such systems 
is split into areas of quasiperiodic motion on phase-space tori like in the integrable 
systems, and into areas where the motion is chaotic. 

The quantum mechanical properties of classically integrable and classically chaotic 
systems vastly differ in the semiclassical limit ?i — > 0. The Wigner functions (the 
quantum mechanical analogues of the classical phase space density) of eigenstates 
of integrable systems are localized on tori in phase space, and the eigenfunctions in 
configuration space have an ordered structure. On the other hand, the eigenfunc- 
tions of fully chaotic/ergodic systems appear randomEP'o* and their Wigner functions 
uniformly cover the whole energy shell in the phase space. 

For the mixed type systems the principle of uniform semiclassical condensation 
(PUSC, see I), I)) states that the Wi gner functions should in this case be localized 
either on the invariant tori in the regions of regular motion or should uniformly cover 
the whole chaotic component of the energy shell in the phase space. This means that 
states are separated into regular and irregular ones. 

In this work we are interested in geometrical and statistical properties of both 
regular and irregular high lying eigenfunctions. For the regular states we employ the 
Einstein-Brillouin-Keller (EBK) quantization, while for the chaotic states we give an 
expression for the wavefunction autocorrelation function based on the geometry of 
the chaotic component in the classical phase space. For more details see referenceEP 
and references therein. Recently we also performed experiments with microwave 
cavities on the same mixed type system as presented hereQ), while theoretical work 
on another mixed type system has been done by Makino et al^ 



* e-mails: gregor. veble@uni-mb. si, robnik(S)uni-mb.si, jliu@tikva.chem.utoronto.ca 
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§2. Definitions 

Our model system was a billiard obtained by conformally mapping the unit circle 
with the complex quadratic polynomial 

z ■w{z) = z + Az^, w{z) = X + iy. 

We chose A = 0.15, where the classical phase space is roughly equally divided into 
components of regular and chaotic motion. The Poincare surface of section (SOS) 
was chosen to lie on the symmetry axis y = with coordinate x and the conjugate 
momentum px as the parameters of the surface. The intersection of the main chaotic 
component of our billiard with the SOS is shown in the figure ||. The coordinate x 
is taken relative to the center of the billiard, while px is the x-component of the unit 
momentum vector. 

The quantum mechanics of billiards is described by the Helmholtz equation 

(Z\ + /c2)^ = 0, (2-2) 

with the Dirichlet boundary conditions, where k"^ = 2mE/h?'. We limited ourselves 
to the states with even parity with respect to reflection across the symmetry line 
y = 0. 

For each state we calculated the smoothed projection of the Wigner function. 
The Wigner function of a state in the general case of N degrees of freedom is 
defined in the full phase space (q, p) as 0* 

^^(q,P) = ^^^ / ^'^Xexp(-ip-X/?i)Vt(q-X/2)V(q + X/2). (2-3) 

In our case the eigenfunctions ip{x, y) generate their Wigner transforms W{x, y,Px,Py) 
through (|2-3| ), where N = 2. In order to compare the Wigner function of a state 
of our system with the classical SOS plot we took its value on the symmetry line 
(y = 0) and integrated it over py, 



Psos{x,Px) = J dpyW{x,y = 0,px,Py). (2-4) 

The result is 

Psos{x,Px) = ^JdX eM-WxX/h)i;^ {x - X/2, y = 0)V(x + X/2, y = 0). 

(2-5) 

Here we see the reason for considering the even parity states only, because ip{x, y = 0) 
is exactly zero for odd states, and therefore a different approach must be used to 
analyze them. 

The catalogue of states studied here consists of 100 consecutive states starting 
at the consecutive index of about 2.5 • 10^. These states were obtained by the scaling 
method first introduced by Vergini and SaracenotiP, that enables us to find a few 
states in the neighbourhood of a chosen wavenumber k. As this is a diagonalizational 
method no levels were missed. Almost each level in our small catalogue can be clearly 
identified as regular or irregular, an idea proposed already by Percivalt 
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§3. Analysis of states 

In a mixed system the phase space is divided into chaotic and regular compo- 
nents. Our work was guided by the principle of uniform semiclassical condensation 
(PUSC, see Robnik 1998), stating that when h tends to the Wigner function of any 
eigenstate uniformly condenses on an invariant object in phase space. This can be 
either a torus in the regular region or a whole chaotic component. Each state could 
thus be labeled as either regular or irregular (chaotic) in the semiclassical limit. By 
looking at the catalogue of states at high (and to some extent even at low) ener- 
gies, one can see that this can indeed be done, though there is still the localization 
phenomenon present due to the still insufficiently low value of the effective Planck's 
constant. 

3.1. Regular states 

We start the analysis by considering the regular states. These are the states 
that can be attributed to quantized tori within the regular regions. For these states 
we tried to employ the EBK torus quantization. We construct a wavefunction on 
the torus as a sum of contributions 



V'i(q) = -4j(q)exp(i 



(3-1) 



of different projections j of the torus onto configuration space. Sj is the classical 
action with respect to some point on the torus and the classical density of trajec- 
tories on this projection. The phase of the wavefunction must change by an integer 
multiple of 2-k when going around any closed contour of the torus. This gives us the 
quantization conditions 



^ / p dq = him + (3i/A) (3-2) 



2tt 

where 7^ are the irreducible closed contours on the torus and Ui the torus quantum 
numbers. The integers [3i are Maslov's corrections and arise due to the changes of 
phase at the singularities of projection of the torus onto configuration space. At 
each caustic encountered along the contour 7j the wavefunction acquires a negative 
phase shift of tt/I, and shifts by vr when refiected from a hard walltl. From this 
consideration it follows that (3i counts the number of caustics plus twice the number 
of hard walls encountered along the contour. 

The task of finding the semiclassical EBK wavefunctions can be divided into two 
parts. The first one is finding the torus with the desired quantum numbers n^, the 
second one being the construction of its appropriate wavefunction in configuration 
space. This is not very straightforward since we do not know the transformation to 
the torus action-angle variables but can only deal witli^the numerically calculated 
orbits. For more details on this procedure see reference . 



* If the contour passes the singularity in the contrary direction to that of the Hamiltonian flow 
on the torus, the phase shifts are of the opposite sign 
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We show an example of a regular state in the figure 2. We present the ex- 
act numerical quantum probability density (top left), the probability density of its 
semiclassical approximation (top right), the classical density of trajectories on the 
appropriate torus (bottom left) and the smoothed projection of the exact Wigner 
function (bottom right). The semiclassical wavefunctions shown are remarkable as 
they possess all of the features of their exact counterparts that are larger than the 
appropriate wavelength. Note that for each torus there are two characteristic wave- 
lengths since there are two quantum numbers associated with it. As it happens in 
our case, the two wavelengths can be of different orders of magnitude. 

3.2. Irregular states 

While for the regular states it was quite straightforward to find their semiclassical 
approximations, the nature of irregular states is very much different. These states 
are very sensitive to small perturbations of the system, so in any physical system 
the individual features of the states are lost when the effective Planck's constant 
tends to 0. The features that are insensitive to small perturbations are, however, 
the statistical properties of spectra and eigenstates. 

One measure of statistical properties of the wavefunctions is the wavefunction 
autocorrelation function, 

r(„ . _ (^t(q^-x/2)V;(q- + x/2))q,,,(q) 

The area of averaging e(q) close to the point q should be taken such that its linear 
size is many wavelengths across, however small enough that the local properties of 
classical mechanics within it are largely uniform. 

If one takes the Fourier transform of the Wigner function ( |2-3| ), it is easy to 
show that 

J W{q, p) exp(ip • x/h)d^p = ^'^(q - x/2)^(q + x/2). (3-4) 

By knowing the Wigner function of an eigenstate, it is then possible to use this result 
to calculate its autocorrelation function. 

According to the principle of uniform semiclassical condensation, the Wigner 
function of any chaotic state should uniformly condense on the whole chaotic com- 
ponent when the effective h tends to 0. Let us limit ourselves only to the cases of the 
Hamiltonians with an isotropic dependence upon p. We can write the semiclassical 
Wigner function as 

Wv, (q, p) = a6{E - H{ci, p))xv, (q, p) (3-5) 

where xVi denotes the characteristic function of the chaotic component and a is the 
normalization constant. 

If we write the characteristic function in two degrees of freedom as a Fourier 
series of the polar angle (j)p of the momentum vector (since the absolute value of p 
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is constant at a given energy and point q), 

oo 

Xx>i(q,p)= (q)exp(im(/)p), (3-6) 



m=— oo 



it is quite straightforward to show by using the above expressions and the integral 
representations of the Bessel functions that 

^ . ^ (Em=-oo'^m (q')^"'^m(p(q')?'/^)exp(im(/):,))q,g,(q) 

CvM.^) = , A/ /x\ ' (^'^^ 

where (px is the polar angle of the vector x. 

We compare the autocorrelation function of a chaotic state with the semiclassical 



prediction 3^ in figure 3. The averaging area e(q) was taken as a circle of radius 0.2 
around the point = (0.65,0). It was taken the same for both the semiclassical 

prediction and for the numerical results. The averaging radius was taken quite large 
in order to reduce the localization properties of the wavefunctions, which are still 
apparent at the values of effective h that we were able to obtain. But this radius still 
has to be taken small enough in order not to completely smooth out the classical 
dynamics. The agreement with the semiclassical prediction is quite good. It clearly 
deviates from the Berry's predicitonffl) for fully ergodic systems and tends towards 
our semiclassical result. 



§4. Conclusion 



In this work we were able to find semiclassical wavefunctions of the regular 
eigenstates, while for the chaotic states we gave a prediction for the wavefunction 
autocorrelation function, both of which match nicely with their exact counterparts. 
What needs to be studied further is the localization of the chaotic states on only 
parts of the whole chaotic component when not yet in the strict semiclassical limit, 
which gives rise to deviations from semiclassical predictions. 
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Fig. 1. The SOS section of the main chaotic component of the A = 0.15 bilhard. 



Fig. 2. The probabihty density for a regular state with = 20421387.1741 (top left), its semi- 
classical approximation with fc^c = 20421385.96 (top right), with the quantum numbers on the 
torus being m = 8648 and n2 = 1 (16 equally spaced contours from to the maximum value 
in both cases). In the bottom row we show the classical density with 20 contours from to 
the maximum of the appropriate torus (left) and the smoothed projection of the exact Wigner 
function (right), with ten contours from to the maximum value. 



Fig. 3. In the top row we show the probability density for the chaotic state with = 20421262.6667 
(left, 8 contours), the circle denoting the region of averaging, with the smoothed projection of 
its Wigner function (right, 10 contours). In the bottom row we plot the wavefunction autocor- 
relation function (averaged as explained in text) in the x (left) and y (right) directions (full), 
compared to the semiclassical prediction (dashed) and the fully ergodic prediction (dotted). 
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